Unraveling hidden hierarchies and dual structures in an integrable field model 
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An integrable field theory, due to path-independence on the space-time plain, should yield together 
with infinite sets of conserved charges also similar dual charges, determining the boundary and defect 
contributions. On the example of the nonlinear Schrodinger equation we unravel hidden hierarchies 
and dual structures and show the complete integrability through a novel Yang-Baxter equation at 
the classical and quantum level with exact solution. 

The Nother's theorem (NT) connects the symmetries of a model to its invariant quantities. In a mechanical system 
with finite degrees of freedom such quantities are the first integrals conserved in time. For the integrability and exact 
£SJ '. solvability of a system the number of conserved charges should match with its degrees of freedom [l| . This definition 
is usually generalized straitforwardly to the field theory (FT) defining the integrability of a field model through an 
infinite set of charges conserved in time. The NT is also tuned to fit this definition, stating that the model should 
have infinite number of symmetries with one to one correspondence to the conserved charges However, in reality 
neither the integrable FT nor the NT insists on such inherently noncovariant properties, which are set perhaps by 
tradition and not by necessity. In fact, the NT relates each symmetry of a field model to an integrated quantity 
invariant along any chosen direction in the space-time manifold Q and similarly an integrable FT, based on the 
flatness condition involving both of its Lax operators, allows a path-independent formulation on the x, t plain with no 
emphasis on the direction of evolution Q. However, as mentioned above, the conventional approach in integrable FT 
focuses mostly on space configurations at a fixed time generating an infinite set of charges conserved in time, through 
the space-Lax operator, which acts as an infinitesimal shift operator (ISO) along the x-direction. This traditional 
formulation however, as we show, works well only for models with periodic or vanishing boundary condition (BC) and 
with no internal defect, but fails for more general BC and in the presence of defect points, where the usual charges no 
longer remain conserved [HH)]. In dealing with the time-dependent boundary conditions, relevant in real experiments 
@, the traditional scattering data along the space direction is found also to be insufficient, prompting the inclusion 
of scattering processes in both space and time directions [4|. These explicit conflicts with the established FT force 
us to rethink about the limitations of the accepted theory and propose a dual approach to complement the existing 
results by unraveling a wealth of hidden structures in the integrable FT. In particular, working on the example of the 
nonlinear Schrodinger (NLS) equation we obtain an infinite set of dual charges conserved in space, generated by the 
time-Lax operator acting as an ISO along the t-direction. This addresses a rather philosophical question, asking how 
£SJ \ many infinite sets of conserved charges one needs for the integrability of a FT, having infinite degrees of freedom. 
By finding an additional set of dual charges ./£ , m = 2, k = 1, 2, . . ., we could enlarge the known set of conserved 

(n) 

charges c k , n = 1, k = 1, 2, . . . for the NLS model by another countable infinity. More encouragingly, we find that, 
the contributions from general BC and the defect point can be given explicitly by these dual charges, which are missed 
in traditional treatments. Extending to higher Lax pairs one could obtain further hidden charges with n = 2,3,... 
along with their duals with m = 1, 3, 4, . . ., generating newer integrable hierarchies covering all possible cases (see Fig. 
1), which are ignored surprizingly in the literature. 

The hardest task however is to show the complete integrability and exact solvability of the system by establishing 
the mutual commutativity of dual charges, in analogy with that for the usual charges, shown by the celebrated Yang- 
Baxter equation (YBE) linked to the space-Lax operator and solved by the well known Bethe ansatz (BA) [8|. For 
this we conjecture a fixed space Poisson bracket (PB) and the related CR for the NLS system, which could derive a 
novel hierarchy of integrable equations from the dual charges and at the same time construct a dual YBE linked to 
the time-Lax operator and solve the model exactly by the BA. 

The integrable NLS equation iq t — q xx + 2(q*q)q, / for the complex field q(x, t) together with its complex conjugate 
can be obtained from the compatibility of linear Lax equations 

$ X (A) - E7i(A)*(A), $t(A) = V 2 (A)$(A), $(A) = (^(A), <fe(A)), (1) 

involving space and time Lax operators U n (X) and V m (X), with n = 1 for the simplest case with space- variable 
x = X\ and m — 2 for the next higher time variable t — £2- In the conventional approach the generating function 
hi(j)i{x, A) = i f^^dx' p(x',\) for the densities p(x,\) = J2k Pk(x)(2\)^ k of charges is considered for vanishing 
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BC using the simplest space-Lax operator Ui(X). Therefore from the first Lax equation in ([I]) one can derive the 
Riccati equation for ^ = T(x, A) = ^ fc T k (x)(2X)~ k , solving which one obtains infinite number of conserved charges 

c k = f^^dx qT k (x), k = 1,2,..., well known for the NLS model Q. The mutual commutativity of the charges 
[cfe, ci] — , fc ^ I is proved by the classical and quantum YBE involving Lax operator ?7i(A) and the rational R- 
matrix, while NLS equation can be obtained from Hamiltonian H = C3 using the usual PB with canonical momentum 
q* Q. The known NLS hierarchy is obtained from higher charges H = c m+1 , compatible with the flatness condition 
of the Lax pairs (Ui, V m ), m = 3, 4, . . . for time t = t m (see Fig. 1). 

For exploring the dual objects we follow a path parallel to the conventional one, replacing however x by t and U\{\) 
by the time-Lax operator V-^A) for the NLS equation. This results from the second Lax equation in Jl}, describing the 
scattering process along the time-direction at fixed x, a generating function ln</>i(i, A) = i J_ ocl dt' j(t' , A) for densities 
j(t,X) — J2k jk(t)(2X)~ k , of the dual charges J k , as j(t,X) — Vu + Vi2T{t,X), involving the time-Lax operator and 
r(t, A) = ~^2 k T k (2X)~ k . Solving the associated time-Riccati equation recursively for we can construct an infinite 
set of dual charges Jk = /_ dt j k (t), fc = 1, 2, . . . , for the NLS equation with vanishing BC on the time-axis (see 
Appendix). Similar idea was used earlier for obtaining nonlocal charges in the sine-Gordon field model Q- 

Importantly, for considering models beyond periodic or vanishing BC in the space interval [xq,Xi], or the models 
with a defect at point xa, the traditional treatment engaging only the space-Lax operator becomes insufficient and 
the inclusion of the time-Lax operator essential. One finds that, the actual conserved charges are not the traditional 
ones mentioned above, but more general ones as 

Ck=c k +I k 3C , Jjf c = J k (x ,t) - J k ( Xl ,t), fc = l,2,..., (2) 

with c k as the traditional charges, while I k given through dual charges J k is the boundary contribution, or similarly 

C k =c-+c++ Jf Jf / = J k (x+) ~ J k (x-), k = 1, 2, . . . (3) 

with vanishing BC but with a defect contribution I^ e ^ given again through J k at = (xa ± e)i e _ > 0) where are 
usual charges for the fields confined to the left (right) from the defect point. It is remarkable that the dual charges 
Jfc, which are derived here but traditionally ignored, determine actually the boundary and defect contributions. 

Switching over to the dual picture we can show similarly that, in more general case of time-interval [0,T] at fixed 
space x, the dual charges conserved in space generalise to 

J k = Jk + T k BC , T k BC = c fc (0) - c fc (T), k = 1, 2, . . . , (4) 
with time-boundary contribution T k c and in models with time-defect at moment t = td to 

Jk = 4 + Jk + T t f > T t f = - Cfc(^), k = 1, 2, . . . , (5) 

with time-defect contribution T k e , where J? are dual charges in the semi-infinite time intervals before and after the 
defect- moment td- 

However, it is natural to ask: how far we can stretch this concept of duality. Can dual charges J k , k = 1, 2, . . . , 
generate a new hierarchy of integrable equations including the NLS equation, mimicking the standard approach 
but using some unusual canonical bracket? Defining the dual canonical momentum as p — j^- we derive from 

the NLS Lagrangian p = q*, p* = q x and conjecture an equal-space canonical bracket {q(x, t),p(x, t )} = 5{t — 
t), {q*(x,t),p*(x,t')} = S(t-t), {q(t),q*(t')} = {p(t),p*(t')} = 0. Using this dual PB and taking J 2 = H as the 
Hamiltonian we derive p* = {p* , J2} = iqt — 2|(?| 2 (7, which by putting p* = q x would yield the same NLS equation 
associated with the Lax pair (C/i, V^)- From next higher order however our dual PB starts producing new equations, 
e.g. for H = J3 one gets a linear equation q t — q X2 related to the pair (U2, V2) and for H = J4, an intriguing integrable 
equation 

iQxt = (lx a + 2(g*g - q*q x ), (6) 

associated with the Lax pair (^3,^2) and so on, which generates a novel integrable hierarchy represented by the 
flatness condition of the infinite series of Lax pairs (U k , V 2 ), k = 1, 2, . . . as a dual to the NLS hierarchy (see Fig. 1). 

Encoraged by the success with our conjecture we set a harder task to establish the complete integrability of the 
system by showing mutual independence of dual charges [J k , J{\ =0, fc ^ I through a novel dual Yang-Baxter relation 
both in classical and quantum cases, where the space-Lax operator U\{X), appearing in the standard YBE for the 
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NLS model [1] is replaced by its time-Lax operator V2(A), while the classical r(A — fi) and the quantum R(X — fx) 
matrices remain the same well known rational matrices (see Appendix). 

However, proving the validity of these dual YBEs becomes exceedingly difficult, due to much more involved structure 
of V^(A) for the NLS model compared to that of Ui(X), . Consequently, in place of just two nontrivial relations appear- 
ing in known YBEs [1,11], ten such equations arise in the dual YBEs, apart from the ultralocal conditions. However 
miraculously the dual PB we propose solves all classical relations exactly, proving the dual YBE for the classical NLS 
model. More significantly, for the quantum NLS model the corresponding equal-space commutation relations (CR) 
can be given for the field operators discretize along the time axis as [q^(x), q* k (x)] — hA^~ 1 Sjk 1 [(p(x), q* (%)] = 0, 
A( t ) being the lattice constant in discrete time. The lattice regularized V 2 J operator together with the rational quan- 
tum i?-matrix and the proposed CR, solve fortunately all nontrivial commutation relations in the dual quantum YBE, 
exactly upto order 0(Am), which is however sufficient for the QFT model of NLS we are interested in (see Appendix). 

Recall that the YBE solutions are the key for exactly solving the eigenvalue problem (EVP) for quantum charges 
Q. Using the same idea we can formulate an algebraic Bethe ansatz (ABA) for the exact EVP solution of the dual 
quantum charges, by constructing a dual monodromy matrix 5(A) = IX/V^CA) with the creation and annihilation 
operators defined through its operator elements as B(X) — 5 i2 (A) and C(A) = 5 2 i(A) respectively, while its diagonal 

elements: t(A) = trS(X) with lnr(A) = Y^ fe Jfe(2A) _fc , are related to dual quantum charges J/., k = 1,2, The 

EVP: r(A)|iV >= A(A)|AT > can be solved exactly in perfect analogy with the standard ABA formulation [8|, since 
the i?-matrix contribution from our dual YBE remains the same. Only the vacuum expectation of U\(X) appearing 
in eigenvalue A(A) in the usual case, should be replaced by the corresponding value for the time-Lax operator V^(A) 
in the dual case. We would like to conclude in general that, the NLS field model posseses rich hidden structures and 
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FIG. 1: Complete NLS hierarchy with Lax pair (U„, V m ), n,m — 1,2, Vertical lines are the hierarchies with infinite charges 

cS-i conserved in time t m and vorizontal lines are are the dual-hierarchies with infinite dual charges J^l\ conserved in space 
x n . The vertical dashed line corresponds to the known NLS hierarchy, while the horizontal solid line to the dual hierarchy 
found here. 



multiple integrable hierarchies beyond the known ones, associated with a double-infinite series of Lax pairs (U n , V m ) 
with flatness condition U ntm — V mXn + [f/„, V m ] — 0, covering all possible cases (n, m) = 1, 2 ... as shown in Fig. 1. 
The diagonal line with n = m represent linear equations since [U n , V n ] = due to U n ~ V n , while the vertical lines 
with m = 1,2... for each fixed n, represent integrable hierarchies with a set of infinite charges c^ +1 conserved in 
time t rn , of which only n = 1 corresponds to the known NLS hierarchy. The hierarchy with n = 2 can be obtained by 
replacing the variables t — > x 2 , x m — > t m in the explicit form of the dual hierarchy found here. The horizontal lines in 
Fig. 1 with n = 1, 2 ... at fixed m, represent on the other hand integrable dual hierarchies with a set of infinite dual 
charges conserved in space x n . m = 1 case corresponds to the equations, where the variables in the known NLS 
hierarchy are replaced as x — ¥ t\, t n — > x n , while the case with m = 2 corresponds to the dual hierarchy found above 
including ([6]). All common crossing points of vertical and horizontal lines must yield the same equations derivable 
from both usual charge and their duals, as we have witnessed above in deriving the NLS equation associated with the 
Lax pair U\, Vi at n = 1, m = 2, from our dual charge as well as from the NLS Hamiltonian . 

(2) 

In summary, using time as dual to the space variable we find explicitly an infinite set of dual charges Jl , k = 1, 2, . . . 

in addition to the set of usual charges in an integrable field model, on the example of the NLS equation. The dual 
charges generated by the time-Lax operator V2 not only determine the defect and the boundary contributions missed 
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in a traditional approach, but also ensure the complete integrability through a dual Yang-Baxter equation, both in 
classical and quantum cases with the use of a novel canonical momentum and an equal-space canonical commutator. 
Similarly, the contributions of time-boundaries and time-defect to the general dual charges are determined in turn by 
the usual set of charges, due to the space-time duality. It is found that, in general Lax pairs (U n . V m ) with all possible 
combinations of (n S N, m 6 N) should generate multiple-infinite sets of integrable hierarchies shown in Fig. 1, only 
one of which (with (n = 1, m G N)) is known, while a few (with (m = 2, n G N), (m G N, n = 2), (m = 1, n € N) ) 
are found by us including the equation ([6]). These are the main results obtained here apart from the exact solution 
of the NLS quantum field model along the dual direction by algebraic Bethe ansatz. The present investigation not 
only exposes the insufficiency in the established integrable field theory, but also complements it and opens up new 
directions of research in integrable classical and quantum systems, like models with time defects and time boundaries 
as dual to such popular problems on the space axis, exploring the NLS multiple-hierarchy found here for their exact 
soliton solutions, Hamiltonian structure, YBE etc [lOj, new class of exactly solvable spin chains along the dual axis 
using different realizations of the quantum time Lax operator, integrable models on an arbitrary path on a 2D plane 
construction of novel quantum algebras as dual to the known Hopf algebras in integrable systems [9] etc. 

Appendix 

I. Infinite set of dual charges 

Using the matrix elements V\\ = 2A + |q| 2 , V\ 2 = V 21 = 2Xq — iq x of the time-Lax operator V2(A) for the NLS 
equation, we can derive its dual charges Jk = J_ dtjk, with jk = {—iQx^k + <7lVi-i)j fc = 1,2,... , through the 
solution Ffc of the time-Riccati equation iY^t — — 2|g| 2 r^ + (gTjTk+i—j — iq x TjTk-j) with densities jk as 

h = i^lq-ffqx), h = i^Q+&x+\q\ 4 , jz = q^Qx, H = «^tfe + 9?9t-«k| 2 (9*ft-9t9)-2|g| 2 q* x qx + (q* 2 ql + qfq 2 ), 
and so on. 

II. Dual classical YBE 

The dual YBE of the NLS model is linked to the time-Lax operator V2(A) with the same rational r-matrix as 
{V(A, x, ti) <g) ! V(/j,, x,t 2 )} — [r(X — fx), V(fx, x, ti) ® I ' + 1 C3) V [n, x, ti)]£(ii — £2), In this matrix relation ten equations 
remain nontrivial, which however using the dual PB can be reduced into only two coinciding or related groups with 
values — iq and — 2i{X + /i) multiplied by <5(ti — t%) on both sides of the equation proving the dual YBE. 

III. Dual quantum YBE 

The dual YBE for the NLS quantum field model with regularized time Lax operator U J and the known rational 
i?-matrix is 

R(X - n)Vi(X) (g> y j (/i) = (I (8 Vi(fj,))(V j (A) ® I)R(X - (x), j = 1, 2, . . .N. 

In the lattice regularization (X) = I + iA( t )V2(A), A( t )| — > 0, the field operators q(x,t),q*(x,t) in the known 
time-Lax operator V2(A) are replaced by their semi-discretize variants q°(x), q* x {x). This results to ten nontrivial 
operator equations containing terms of order 0(A.u\) apart from higher order terms. For example, in element (11, 12): 
V(VAWV? 2 (/;) - V( x {p)Vl 2 (A)) + (A - /i) (Vl {X)V( 2 ( M ) - V( 2 (fi)V^ (A) ) , the first term is nontrivial: -2r?A (t) (A - fit? . 
However the second term containing i A 2 t ^ ( A — pt) [qi \ 3 1 , q x ] , cancels it exactly by using the proposed dual CR. Similar 
picture is repeated for other terms inducing exact cancellation of all terms upto order O(A^). The ultralocal condition 
VHX) ® V k {fi) = V k (X) ® 7)(7 ® F^'(/i)), k + I. follows from the canonicity of the dual CR. 
The author is indebted to Prof. P. Mitra for stimulating discussions. 
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